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Abstract 

This paper deals with several technical issues of non-perturbative four-dimensional 
■ Lorentzian canonical quantum gravity in the continuum that arose in connection 

with the recently constructed Wheeler-DeWitt quantum constraint operator. 

1) The Wheeler-DeWitt constraint mixes the previously discussed diffeomorphism 
q ■ superselection sectors which thus become spurious, 

2) Thus, the inner product for diffeomorphism invariant states can be fixed by re- 
quiring that diffeomorphism group averaging is a partial isometry, 

3) The established non-anomalous constraint algebra is clarified by computing com- 
mutators of duals of constraint operators, 

4) The full classical constraint algebra is faithfully implemented on the diffeomor- 
phism invariant Hilbert space in an appropriate sense, 

5) The Hilbert space of diffeomorphism invariant states can be made separable if a 
natural new superselection principle is satisfied, 

6) We propose a natural physical scalar product for quantum general relativity by 
extending the group average approach to the case of non-self-adjoint constraint op- 
erators like the Wheeler-DeWitt constraint and 

7) Equipped with this inner product, the construction of physical observables is 
straightforward. 
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1 Introduction 

Recently |], 0, [|, a new technique was successfully employed to quantize the Wheeler- 
DeWitt quantum constraint operator of four- dimensional Lorentzian canonical gravity in 
the continuum (see also H] for a nice and compact non-technical exposition of the essential 
results of these papers). More precisely, the results of these papers show that while using 
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all the rigorous kinematical framework that has been developed in Jj[ || [7|, 0, |10j, it is 
actually possible to define a densely defined operator on the underlying kinematical Hilbert 
space whose classical limit corresponds to the classical Wheeler-DeWitt (or Hamiltonian) 
constraint of canonical general relativity. 

The Hilbert space itself provides a representation of the operator algebra in which 
an SU(2) connection is diagonal, that is, it uses very crucially the fact that we are 
dealing with a connection dynamics rather than geometrodynamics formulation of general 
relativity. The advantages of the use of a connection as opposed to a three-metric has 
been discovered a decade ago by Ashtekar However, in contrast to the original point 



of view advocated in [11] and many later publications which made essential use of a 
complex valued connection, we are dealing with an entirely real formulation of Lorentzian 
general relativity, a point of view first stressed by Barbero ]12[ . Although the Hamiltonian 



constraint of general relativity reads much more complicated in terms of real than in terms 
of complex connections, it allows us to solve the difficult reality conditions which arise 
for the complex variables and which were one of the two major roadblocks to making 
progress with these (see, however, |TS|] for partial progress in this direction). Moreover, 
and most importantly, one is able to rigorously quantize the loop variables introduced by 
Gambini and Trias and rediscovered for general relativity by Rovelli and Smolin |Tj| , 



by means of techniques developed in 0, || [7|, [|, [9|, |10fl . 

The virtue of |, g, 1 is that one can overcome the second major roadblock to quantiz- 
ing the Wheeler-DeWitt constraint : its highly non-polynomial structure in terms of the 
basic variables which make seem it impossible to give corresponding quantum analogs any 
mathematical meaning (one can make the Wheeler-DeWitt constraint polynomial upon 
multiplying it by a density weighted overall factor but this is not allowed in a diffeomor- 
phism invariant theory for reasons explained in [I], and in fact leads to a highly singular 
operator). The results of |L], 0, |] also show that the operator so constructed is anomaly- 
free in a suitable sense and in [Q the entire space of rigorous solutions to all constraints 
was constructed. However, the quantization programme was still not completed because 
the following questions were still unresolved : 



1) In ||T0| the existence of superselection sectors for diffeomorphism invariant theories of 
connections was dicovered (theories which do not have any further constraint). The 
superselection rule holds provided that one can deal with strongly diffeomorphism 
invariant observables only. Is this assumption justified in general relativity ? 



2) In 1X0 1 an infinite number of different inner products, appropriate for the general 



solution to the diffeomorphism constraint, all of which implement the reality con- 
ditions, were constructed. An infinite number of overall factors were kept undeter- 
mined for the inner product on each sector. Does the presence of the Hamiltonian 
constraint somehow fix those factors ? In other words, is there a unique group 
averaging procedure |23|, |2^] (employed there to construct an inner product on dif- 
feomorphism invariant states) ? 

3) The operator constructed in [|I|, [2| suffers from a huge amount of ambiguity : There 
is not only one operator but an uncountably infinite number of them, the freedom 
being captured by a regularization choice. The dual of that operator, determined by 
evaluating it on a diffeomorphism invariant distribution, however, was shown to be 
choice independent because the choice freedom is equivalent to a choice of (partial) 
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diffeomorphism gauge fixing. The problem then, however, was that the dual of 
the Hamiltonian constraint is not diffeomorphism invariant and so multiplying dual 
operators was not obviously defined. Can this be clarified ? 

4) In [JJ] the commutator of two Hamiltonian constraint operators was shown to vanish 
when evaluated on diffeomorphism invariant states. However, we did not compute 
the commutator directly between duals of these operators. More seriously, it is 
unexpected that the commutator just equals zero instead of being proportional 
to a diffeomorphism constraint operator. Finally, we also did not compute the 
commutator between the Hamiltonian and the Diffeomorphism constraint operators. 
Can these problems be at least partially solved ? 

5) The Hilbert space of diffeomorphism invariant states is non-separable in the present 
framework. Can one make it separable by imposing a supersection principle ? 

6) In || an algorithm was derived of how to find the general solution to all quantum 
constraints in principle (to find such solutions is a matter of computations and can 
be done in finite time). However, we did not provide an inner product on the space of 
solutions which would separate a general distribution from a normalizable physical 
state. The main problem is that the Hamiltonian constraint is, expectedly, not a 
self-adjoint operator and so the group averaging technique which was designed for 
self-adjoint constraint operators cannot be applied. Is it possible to supplement the 
solution space with an inner product ? 

7) Since in || no physical inner product was derived, the issue of Dirac observables 
was unresolved. Can we construct such observables ? 

The aim of this paper is to answer all these questions. Specifically we show that : 

1) The superselection sectors are spurious : the Hamiltonian constraint mixes the sec- 
tors and therefore solutions are linear combinations of vectors from different sectors. 
Interesting physical observables therefore cannot leave these sectors invariant. 

2) Since the superselection sectors are spurious the group averaging maps on each sector 
must come from one big averaging map. We show that there is only one such map 
if one wants to implement the group averaging such that strongly diffeomorphism 
invariant observables are self-adjoint and such that a natural class of orthonormal 
bases on the kinematical Hilbert space stays orthonormal after averaging. 

3) As was shown in ||, arbitrary products of Hamiltonian operators have a dual which 
is again choice independent when evaluated on dioffeomorphism invariant states. 
Thus, one can define the product of duals to be the dual of the product in reverse 
order. Therefore commutators can be rigorously computed and trivially reproduce 
the results from @]. 

4) As foolows immediately from |1 the commutator algebra between dual Hamiltonian 
constraint operators is Abelian which is not expected. However, we prove that there 
exists an operator 0(M, N) corresponding to / d 3 x(MN^ a — M^ a N)q ab Vt, (which is 
the right hand side of the Poisson bracket {H(M), H(N)}) and whose dual is again 
choice independent and annihilates diffeomorphism invariant states. We therefore 
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have the trivial operator identity [H'(M), H'(N)] = —0'(M,N) on diffeomorphism 
invariant states (here the prime stands for dual) and it is in this sense that the Dirac 
algebra is faithfully implemented. Morover, the commutator between Hamiltonian 
and diffeomorphism constraint closes in the expected way. 

5) We argue that the present results and the structure of the Hamiltonian constraint 
(dynamics) suggest that the Hilbert space breaks into a continuous orthogonal sum 
of isomorphic sectors each of which is labelled by a countable number of moduli 
parameters and that no physical observable maps between those sectors. Thus, 
physically one may restrict to one of these sectors which turns out to be separable. 

6) We extract the essential features of the group average map for self-adjoint operators 
and show that these can be naturally extended to any non-self-adjoint operator. 
We show that in quantum mechanical cases our proposal does lead to the known or 
expected result. Astonishingly we find that even for self-adjoint operators the group 
average proposal can give rise to solutions (obtained via the rigorously defined group 
averaging map) to the constraint which are in fact elements of the original Hilbert 
space, in spite of the fact that the constraint operator is unbounded ! We provide 
quantum mechanics models which mimic the Euclidean or Lorentzian Hamiltonian 
constraint and find that all solutions are normalizable with respect to the original 
kinematical inner product. This suggests that there is nothing wrong with the fact 
that the basic solutions found in || are in fact all normalizable with respect to the 
inner product obtained by group averaging the diffeomorphism constraint. 
Finally, we use our method to propose a physical inner product for quantum general 
relativity which is very natural. It is very different from the kinematical diffeomor- 
phism invariant inner product and encodes the full dynamics of quantum general 
relativity. 

7) We take the definition that a Dirac observable maps the physical Hilbert space to 
itself and is self-adjoint thereon. Given the inner product constructed, the construc- 
tion of a complete set of at least symmetric operators is straightforward. 

We will leave unresolved the very physical questions of quantum general relativity and 
leave it for future publications. An (incomplete) list of such questions is given by : 



1) The problem of time, deparametrization of the theory (see [16] for first attempts in 
this direction). 

2) Physical interpretation of Dirac observables or how to represent known classical ob- 
servables in the present framework. 

3) Existence of a semi-classical limit |25 . 



4) Black holes and the entropy puzzle from first principles (see [|18| , |17| , |T9|j for first success 
in this direction). 

5) Connection or relation with string theory. 

6) Matter coupling and Energy operator J2(], 

The organization of the article is that we answer the questions raised above in the chrono- 
logical order displayed after introducing some basic notation at the beginning of the next 
section. 
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2 Diffeomorphism Superselection, Graph Symmetries 
and Diffeomorphism Invariant Scalar Product 

The goal of this section is a) to show that there are no superselection sectors in quantum 
gravity once we take the Hamiltonian constraint into account and b) to select an in- 
ner product for diffeomorphism invariant states of quantum gravity. On the other hand, 
the superselection sectors arose by restricting the algebra of observables to those which 
strongly commute with all diffeomorphisms. We show that if we require that these observ- 
ables are to be promoted to self-adjoint operators and that group averaging with respect 
to diffeomorphisms is a partial isometry, then there is a unique diffeomorphism invariant 
group averaging inner product. 

We begin with a compact review of the relevant notions from [[TIJ. The interested 
reader is urged to consult this paper and references therein, in particular H [| [7], [| . 

By 7 we will denote in the sequel a closed, piecewise analytic graph embedded into a 
d-dimensional smooth manifold £ (the case of interest in general relativity is d — 3). The 
set of its edges will be denoted E(j) and the set of its vertices ^(7). By suitably subdi- 
viding edges into two halves we can assume that all of them are outgoing from a vertex 
(the remaining endpoint of the so divided edges are not vertices of the graph because 
they are points of analyticity). Let A be a G-connection for a compact gauge group G 
(the case of interest in general relativity is G = SU(2)). We will denote by h e (A) the 
holonomy of A along the edge e. Let itj be the (once and for all fixed representant of 
the equivalence class of the set of independent) j-th irreducible representations of G (in 
general relativity j is just a spin quantum number) and label each edge e of 7 with a label 
j e . Let v be an n-valent vertex of 7 and let ei, .., e n be the edges incident at v. Consider 
the decomposition of the tensor product ®~k = i^j £k into irreducibles and denote by c v the 
linearly independent projectors onto the singlets that appear. 

Definition 2.1 A spin-network state is defined by 

T 7 J,s( A ) '■= tr(® ve v(r,)[cv ■ ®eeE( 7 ),veeKj e (he(A)))) (2.1) 

where j = {j e }e£E(j), c = {c v } ve yM. In what follows we will use a compound label 
7^(7(1), /(/),?(/)). 

Thus, a spin-network state is a particular gauge invariant function of smooth connections 
restricted to a graph. Their importance is that they are an orthonormal basis for a Hilbert 
space 7i = Ti. aux , called the auxiliary Hilbert space. Orthonormality means that 

T J gi 2y 7/ j» ^ > aux = Oyy'O'jjiOcjff ■ 

Another way to describe TC is by displaying it as a space of square integrable functions 
L2{A/Q,dfio). Here A/Q is a space of distributional connections modulo gauge transfor- 
mations, typically non-smooth and fio is a rigorously defined, a-additive, diffeomorphism 
invariant probability measure on A/Q, called the Ashtekar-Lewandowski measure. The 
space A/Q is the maximal extension of the space A/Q of smooth connections modulo 
gauge transformations such that (Gel'fand transforms of) spin-network functions are still 
continuous. The inner product can be extended, with the same orthonormality relations, 
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to any smooth graph with a finite number of edges. It can actually be extended to the full 
smooth category |26j but only at the expense of introducing a huge amount of additional 
technicalities. 

We will denote by <3? the finite linear combinations of spin-network functions and call it 
the space of cylindrical functions. A function / 7 is said to be cylindrical with respect to a 
graph 7 whenever it is a finite linear combination of spin-network functions on that graph 
such that TTj e is not the trivial representation for no e G E(j). The space $ is equipped 
with some topology, for instance the folowing "Fourier topology" ||/||i = J2i I < Tj,f > | 
(called this way for reasons explained in ||27|| ) which turns it into a topological vector 
space. By $' we mean the topological dual of $, that is, the bounded linear functionals 
on $. By the Schwarz inequality we have the inclusion $ C K C (notice, howevver, 
that this is not strictly speaking a Rigged Hilbert Space because the Fourier Topology is 
not nuclear, see |10| for another definition). 

This will be enough background for the purpose of the present paper. 



We now begin entering the diffeomorphism invariant regime. 

Definition 2.2 i) The set of smooth diffeomorphisms which leave a piecewise analytic 
graph 7 piecewise analytic is denoted by Diff^CE). The set of all smooth diffeomorphisms 
will be denoted Diff'(H). 

ii) Let fry be a function cylindrical with respect to a graph 7 and <p G Diff(Y,). A unitary 
representation of Diff(Y?) is defined by U(ip)f^ = f^)- We define also U^Tj =: T v .j. 
Hi) A srongly diffeomorphism invariant operator O onTi is a self-adjoint operator on H. 
such that U(ip)0 = OU(ip) for any if G Diff{T). 

Notice that in contrast to [|H| we do not work with analytical but with smooth diffeo- 
morphisms. 

The naive idea of how to obtain a diffeomorphism invariant state is to take an element of 
$ and to consider its orbit under diffeomorphisms. More precisely we have the following 

Given a cylindrical function / G $ we can uniquely decompose it as / = fji where 
/ 7i is a function cylindrical with respect to the graph 7^ and all the graphs in the sum are 
considered to be different. The first guess is to consider the orbit {/} of the state / under 
diffeomorphisms, that is, all cylindrical functions /' such that there is a diffeomorphism 
ip with U((p)f = f and then to just define [/] := J2fe{f} /'■ We immediately see that 
this does not work : 

To see this consider the simple state / = Tj + Tj where 7(f), 7( J) are disjoint, say. Then 
there are infinitely many analyticity preserving diffeomorphisms ipj to be considered in 
the orbit {/} which leave 7(f) invariant but not 7 (J) and vice versa. As a result, [/] 
would contain the meaningless term l]Tj. The reason for this difficulty is that each 
graph has an infinite-dimensional invariance group of diffeomorphisms. Therefore, func- 
tions cylindrical with respect to non-coinciding graphs have to be averaged separately. 
Thus we are lead to define 

N N 

[/]:=E E /[ = £[/»] 

where clearly [/ 7 ] only depends on [7], the orbit (or generalized knot class) of the graph 
7- 
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However, this still does not quite solve the problem because the orbit size of different 
states cylindrical with respect to the same graph can be different. To see what the prob- 
lem is, consider the following instructive example^] : Let 7 be the figure-eight loop (with 
intersection) and consider two spin-network states on 7, T\ := T 7i j lj;2>c with ji 7^ j 2 and 
T 2 := T 7i j JiC / where c, d are some contractors and ji, J2, J are certain irreducible represen- 
tations. Now notice that there is a diffeomporphism ipo (the one that rotates the graph 
by 180 degrees) which leaves T 2 invariant but not T\. As a result we see that the orbit 
size of Ti is double as large as that of T 2 . Therefore we find that [Ti + T 2 ] = [T x ] + 2[T 2 ] 
so that the [.] operation is not at all linear ! 

Since we want a linear operation as otherwise we will not preserve the linear structure 
of quantum theory, we must modify the definition of [/ 7 ] even further, namely, we must 
choose a basis, write f 1 in terms of this basis and average each basis element separately. 
Thus, the group averaging becomes basis dependent, for each basis we get another average 
map. 

Not every basis is allowed : we want the resulting average map to commute with the dif- 
feomorphism group. Notice that this is not an entirely trivial requirement because given 
a basis element Bj and a diffeomorphism (p we are supposed to decompose U((p)Bj into 
basis elements Bp before averaging and then to average the Bp. Thus, it is not clear 
that [U((p)Bj] is the same as U(<p)[Bj] = [Bp]. We will call a basis with this property an 
allowed basis. 

This basis dependence is an ugly feature because it implies that unitary operators are 
not promoted to unitary operators. However, we will see that the Hilbert space splits into 
two orthogonal subspaces Ti s ,Ti A which are basis independently defined : Namely, Ti A 
is nothing else than the kernel of every group averaging map while Ti is its orthogonal 
complement which shows that both spaces are invariantly defined. In each of the subspaces 
the group averaging can be defined basis-independently while being a linear a map. In 
other words, there are preferred allowed bases of Ti which respect the orthogonal split 
Ti = Ti s © Ti A and the group average with respect to either of them is equivalent. As 
one might expect from the above analysis, the spin-network basis is not a preferred basis, 
although it is closely related to one. 
Let us make precise what we have just said. 

Definition 2.3 i) Let {Bj}j be an orthonormal basis for Ti (for instance a spin-network 
state basis with Bj = Tj, without loss of generality we can assume that still I = (7(1), A(/)j. 
Consider the orbit of a basis element Bj under Diff that is, 

{B,} := {U{ V )Bp V e Diff l{I) {Y)}. 
The group average of Bj is defined to be 

B m :=[B x ]:= £ B G <£' (2.2) 

Be{Bj} 

and the group average of any f = Y.k=i a kBi k G $ with respect to the basis {B IyI is 
defined by [f] := EfcLl a fc-%J- 

ii) A basis {Bj} is said to be allowed if a) the corresponding group averaging commutes 
1 This example is due to Don Marolf. 
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with the diffeomorphism group, i.e. [U(ip)Bi] = U((p)[Bj] = [Bj] for all I,ip G Diff^rjJTi) 
and b) all the vectors are well-defined distributions, that is, [Bj\ G $' . In the sequel we 
will only allow for average maps determined by allowed bases and without loss of generality 
we restrict to orthonormal bases. 

Hi) The diffeomorphism invariant inner product with respect to group averaging as deter- 
mined by an allowed basis {B{\ is defined to be 

< [Bj], [Bj] >Diff.= [Bi\{Bj) := J2 <B, Bj > aux (2.3) 

Be{B t } 

and extended by sesquilinearity. We will call the set of group averaged cylindrical functions 



Notice that in contrast to [10| we do not average with respect to the graph 7(f) but with 



respect to the whole state labelled by /. This difference will have important consequences 
with respect to the so-called graph-symmetries to which we turn in a moment. Also notice 
that we average with respect to analyticity preserving diffeomorphisms rather than with 
respect to analytic ones : while this reduces the size of diffeomorphism invariant states (it 
increases the orbit and two analytic knot classes may now lie in the same orbit) it poses 
no problem and all the properties derived in Jl0| are preserved with the bonus that we 



do not have to deal with the subtle type I and II graphs any longer that were discussed 
there. 

We now construct the subspaces H, s ,li, A mentioned above. 

Definition 2.4 i) Given a graph 7 the set of graph symmetries of '7 is the subset G(j) C 
-Dz/f 7 (X) of diffeomorphisms which leave 7 as a whole invariant but permute its edges and 
vertices. 

ii) Given the orbit {Tj} of a spin-network state Tj choose a set D({Ti}) of diffeomorphisms 
which are in one-to-one correspondence with the orbit. Let ipi, ip 2 G D({Tj}) and say that 
they are equivalent iff there is a ip G G((p± (7OO)) such that <p 2 = V^o fx- We define 
D ({Tj}) to be the set of equivalence classes of these diffeomorphisms. Clearly D ({Ti}) 
is in bisection with D([-f(I)}) where [7] is the orbit of the graph 7 under Diff^(E) and 
-DQ7]) is a set of diffeomorphisms in bisection with [7]. 

Hi) Let Tj be a spin-network state. Let us write I = (7(1), A(/)). Then define 

A([/]) := {A, 3<p G G( 7 (/)) 3 U(<p)Ti = T 7(7) , A } and n([I}) := \A([I})\ . (2.4) 

To see that A [J] interacts well with the diffeomorphism group we need the following 
lemma. 

Lemma 2.1 The set A([J]) is diffeomorphism invariant. 
Proof : 

Notice that if ^0(7) = 7 then ipoip o(p~ l ((p(p{)) = ^(7) for any ip G Diff 7 (S). We conclude 
that G(ip(j)) = ip o G(j) o {p^ 1 . Let ip' G -D({T/}). Then we may pick ip G -D ({^/}) 
such that U(ip)Tj = T^^j)^x(i)- Since there is <p G G(ip'(j(I))) = G(<^(7(/))) such that 
ip' = ipo o (p we see that A ([<£>' ■ I]) = A ([/]). 
□ 

The existence of a non-trivial graph-symmetry group implies that certain different la- 
bellings of the graph are related by diffeomorphisms. 
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Example : 

Consider, again, the figure-8 loop 7 = a\ o a<i- Label the loop ctj with ji, % = 1,2 and 
the vertex with c. The figure-8 loop has a graph symmetry-group consisting of <£>o, id 
where <p (ai) = ay, i ^ j and thus A([i]) = {(j'i,j2,c), (j 2 , ji,c)} provided j 1 ^ j 2 and 
A([/]) = {(j, j,c f )} if h=h=j. 

We see that the orbit size of the state depends not only on the graph but also on the 
labels. It is for this reason that the spin-network states are not a well-adapted basis for 
the group-averaging with respect to diffeomorphisms. For example, consider a strongly 
diffeomorphism invariant operator O and the states ipi — ^yjij^c^ = ^twji.cj 03 = 
T i,j,j,c with Ji h as above with U{(p )ipi = tp 2 , ^(^0)^2 = 0i, £%o)^3 = fa- 
First of all we see that [ipi] = [^2] (here we are tentatively assuming that the spin- 
network basis is an allowed basis and average with respect to it). Moreover, suppose that 
0*01 = aipi + bip 2 + eip 3 , 
Oijj 2 = bvpi + aip 2 + c0 3 , 

O0> 3 = dipi + dip2 + eip3 with arbitrary complex coefficients a, b, c, d, e. 
This is the most general strongly diffeomorphism invariant operator on these three states. 
If it is self-adjoint then a, b, e are real and c = d. Now let us see whether with the present 
definition O is promoted to a self-adjoint operator on Ti^ff. The necessary and sufficient 
condition is that [ipi\{Oipj) = [O^]^-). But [ipi](Oip 3 ) = 2d because [ipi] = [1^2] so that 
the orbit of also contains ip2 while \Oip\\{ip^) = c = d contradicting symmetry. 
On the other hand, this example shows how to cure the situation by a suitable change 
of basis : consider the orthonormal states ip± := ± ^2) /V%- Then by definition 
[ifj + ] = v^fV'i], [^-] = 0, that is, the group average as defined by (|2.2j ) has a huge kernel 
even when restricted to vectors defined on the same graph. Notice that group averaged 
states remain orthonormal if not zero, provided we divide the average of ip + by the number 
of spin-network states involved (here two) since they come out of the averaging anyway : 
|[V ; +](V'+) = ^gh^iK^i + 1P2) = 1 (still averaging with respect to the spin-network basis) 
and the so defined inner product is hermitean. In terms of the new vectors the operator 
O becomes 

00+ = ^+ + >/2c0s, 
00- = ffl-, 

Oip3 = d^/2ip + + eipz and we see that it leaves the kernel of the group averaging map 
invariant and moreover does not have any matrix elements between the kernel and its 
complement. We now verify that ^[ip + ](Oip 3 ) = -4= = -A= = ~ [Oip + ] (ip 3 ) is indeed self- 
adjoint. How did this happen ? It happened because the group average was defined 
with respect to the spin-network basis and so one must expect a basis dependence of 
the resulting diffeomorphism invariant inner product. Indeed, the average with respect 
to the new basis ip± modified the averaging as compared to the spin-network basis and 
in fact, the spin-network basis is not well- adapted to difeomorphism invariance because 
the group average of two orthogonal states with the same graph on TC can result in 
diffeomorphism invariant states which are identical and so have unit inner product. By 
changing to the basis involving ip± we have done two things : first, the symmetric state 
ip + is essentially what is important after group averaging while the anti-symmetric state 
■0__ is unimportant and becomes zero after averaging. Morover, the representation of 
diffeomorphism invariant self-adjoint operators is completely reducible with symmetric 
and anti-symmetric states as invariant subspaces. 
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We will now show that these features of the example are true in general and that the 
generalization of the basis of the ip± is preferred. 



Definition 2.5 i) Given a spin-network state Tj, label the n = n([I]) states T 7 m t x k '■= 
T/ ;fc where {X 1 , .., \ n } = A([J]). Now construct the orthonormal states 



1 k 

T t ■= /,„ lE^- fcT Wi]Mfc=l,,n-l. (2.5) 

yjk(k + l) l=l 

ii) Let 7i s be (the completion of the) span of the Tf and likewise for H A . Then 7i = 
HP © H A . 7i A is called the kernel of the average map [for reasons explained in the next 
theorem]. 

Hi) A basis such that each basis vector lies either in H s or in H A is said to be adapted. 

That these states are orthonormal is readily verified by induction over k = 1, .., n. We see 
that we arrive at a new orthonormal decomposition of the Hilbert space of states and of 
the space of functions cylindrical with respect to given graph 7. 

Theorem 2.1 i) The orbit size of any state of the form Tf is given by D([y(I)]), that 
is, it only depends on the orbit size of the graph and not on the state, 
ii) The adapted basis given by the Tf and T A k is allowed. 

Hi) The space [7i A ] 7 is the kernel of every average map restricted to functions cylindrical 
with respect to a graph 7 for any 7 and therefore is invariantly defined, 
iv) Two average maps defined by adapted bases ofTi are identical. We will call the resulting 
average map rjoiff- 

Proof : 
i) 

since the Tf are invariant under the graph symmetry group G("y(I)) they actually do not 
depend on X(I) but only on A([J]). Therfore, the orbit size of Tf is that of the set D ({T/}) 



which in turn is in bijection with D([y(I)]). Notice that therefore rjDiffTf = [T/,i]/ yn([I]) 
where [T/^] is the average with respect to the spin-network baisis. This shows that 

(vDi ff Tf)(Tf) = 1 [r J>1 ](r J , 1 + .. + r J>B([J]) ) 

^n([I\)n([J\) 

n([J}) lT ]( , _ n([J]) 



^n([I])n([J]) ' " ^n([I])n([J]) 
=■■ S im (2.6) 

so the states rjoiffTf are orthonormalized under the associated diffeomorphism invariant 
inner product. 

ii) 

We need to show that [Tf] = \U(^)Tf} and [T A k \ = = [U(<f)T A k ] for each <p, I, k. How- 
ever, by definition, each of the Tf, T A k is nothing else than a particular linear combination 
of spin- network states Tj which have the nice feature that a diffeomorphism just can be 
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translated into an action on the label / : U((p)Tj = T^^^i) = T^-i- It follows that 

U(ip)Tf = Tfj and U{^)Tj A k = T A T ^. k . Thus the result follows from the fact that Tf 

only depends on A(J) which is diffeomorphism invariant and from the fact that [T A k ] = 

for any /, k. 

iii) 

=>• 

First we show that every state of the form Tj — U((p)Tj, where cp £ G(7(J)), is averaged 
to zero by every averaging map, that is, by the averaging map determined by any allowed 
basis {Bj}. Because the averaging map is linear and these states span TC A , we will have 
shown that Ti A is in the kernel of every average map on TC A . 

Since Tj and Bi belong orthonormal bases there is a unitary operator V with matrix 
elements Vjj such that Tj = J2j VijBj and the sum extends only over labels such that 
7(J) = 7(/). Now the state U(tp)Bj is a linear combination of states Bj> : in other 
words U(tp)Bj = J2j> Ljj'Bjr. Therefore by definition [U(tp)Bj] := Ljj>[Bji] = [Bj] 
because {Bj}i is an allowed basis. We conclude that [Tf] = [U(if)Tj]. 
<= 

Conversely, we need to show that every vector that is averaged to zero and that is cylin- 
drical with respect to the same graph 7 lies in Ti A . Since we have shown already that each 
vector in Ti. A has this property we need to show that no vector in has this property 
because the average map is linear. 

Let us introduce a new notation : Let Tf := Tf and T A k := T A and consider any allowed 
orthonormal basis Bj. There exists a unitary operator V such that 

Tf = ^ VuBj, T A = V^jBj and Bj = £ %{Tf + £ V,jT A . 

II i H 

According to i) all vectors Tf are invariant under the graph symmetry group (7(7) while 
the vectors T A are not. In fact, by definition of the graph symmetry group, their orbit size 
is an integral multiple of the orbit size of 7 which is the orbit size of any Tf. Therefore the 
orbit size of the state Bj is an integral multiple n(I) > 1 of the orbit size of 7 the value 
of which depends on which and how many of the states T A occur in its decomposition. 
Combining this with the first half of the proof of iii) we find \Bj\ = n(I) J2i VnVDiffTf . 
Certainly, n(I) is a bounded function of / as otherwise not all the [Bj] would be in 
(since [vDiffTf](Tf) = 5 A{[im[J]) ). 

We now ask whether there exist complex numbers Oj such that the vector ip s := J2i a iTf £ 
7i is averaged to zero. By definition we have 

i I i I j 

= E VmjfTf E OiE vMi)v jT ] = 

j i I 

which has a non-trivial solution {cii}i if and only if the operator M on 7Y' 7 with matrix 
elements between the Tf 

M lJ = Y.Vun(I)V 3l =[VbV% 
1 

is singular where we have defined the operator t) on 7i defined by DBj = n(I)Bj, that 
is, it has matrix elements D LJ = n(I)5 L j between B T , Bj. Notice that the operators D, V 
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are bounded so that the resummation that we performed above is justified. 
Since n{I) > 1, the operator D has a square root R and thus M = [VR][VR]^. Thus we 
are asking whether the projection to 7i s of a positive definite operator on TC is singular 
which is impossible. 

This shows that [ip s ] = if and only if ip s = 0. 
iv) 

Let {Bj}j be an allowed basis adapted to the orthogonal decomposition TC S © H A . That 
is, each Bj lies entirely only in one of these spaces. Since the Bj define an orthonormal 
basis there is a unitary map V" 5 © V A such that (V S )W A = {V A )W S = 0, {V S ) ] V S = 
ids, (V A yV A = id.A and Tf = J2j VfjBj, T A = J2j V^jBj where the sum runs only over 
labels J such that 7(f) = 7 (J) (we have written T A k = T A ) and also only Bj are involved 
which are either in 7i s or 7i A respectively. By definition 

P/] = EVS E Xj^vu E u{<p)Bj 

J B'jE{Bj} J tp<=D({Bj}) 

where, as before, D({Bj}) is a set of diffeomorphisms which are in one to one correspon- 
dence with the orbit of Bj. Notice that Bj = Ek(V S )jk t k 

since the basis is adapted 

(there are no T A on the right hand side). According to i) the orbit size of all the Tj| is 
identical and in fact coincides with the orbit size of the graph 7(1). That means that 
each if G D({Bj}) which moves Bj at all must move all the T% which appear on the right 
hand side of Bj = Y,k(V S )jkTk- ft follows that the orbit size of Bj is the same as that 
of [7(1)] as well. It follows that [Bj] = Y J K{^ S )~jK r lDif fT§ so that, putting these results 
together, we find 

[T/] = nmffJf 

where r]mff is the average map defined by the basis Tf . 
□ 

Let us verify that r) Di ffTf E If / G <&, then 

\[Ti S ](f)\ := -7^=1 E <TJ> aux \< n L= \<TJ> aux \ 
ln(I) Te{T 7 } \M J ) ^mi 



< -7^=El< T ^/> a -l = 4== ( 2 - 7 ) 



which shows that r] Di ffT b G ^'m/f- 

Let us summarize the contents of Theorem ( |2.1|) with regard to the uniqueness of the 
average map : 

Every average map has the same kernel Ti A so that the kernel is defined independently of 
the allowed basis that was used to define the average map. Therefore the spaces TC S , TC A 
are invariantly defined. 

Being invariantly defined one can consider all bases that are adapted to the associated 
orthogonal decomposition which is natural because the space 7i A does not play any role 
after moding out with respect to the diffeomorphism group. Moreover, as we will see be- 
low, the subalgebra of the observable algebra given by strongly diffeomorphism invariant 
observables leaves the subspaces TC S , TC A invariant so that we have an additional motiva- 
tion to consider only averages with respect to adapted bases in the sequel. 
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Then, when restricting to adapted bases it turns out that there is actually only one aver- 
aging map. 

Thus, among all the possible averagings, one of them, ^Di//, is singled out. 

We could be content with this result, however, we now ask if we can do even better : 
Recall that we defined an average map by first decomposing a cylindrical function into 
components each of which was cylindrical with respect to a certain graph. Now diffeo- 
morphic graphs certainly have identical orbit sizes but non-diffeomorphic ones do not 
have, in general, the same orbit size and as we have seen above, in general their orbit 
sizes are "infinitely different" in the sense that one gets in general not an element of $' if 
summing over all vectors in the orbit of a vector which is a linear combination of vectors 
cylindrical with respect to non-diffeomorphic graphs. It seems therefore that we do not 
have any justification for a universal average map rjmff, that is, for each diffeomorphism 
equivalence class of graphs [7] we could have a different average map Vmff- Acoording to 

the above analysis, the difference between rjoiff and /To!// can really be only a pre-factor 
if we choose the canonical way of averaging and the prefactor is a positive real number 
which, heuristically speaking, accounts for the fact that each graph has different orbit 
size as compared to the orbit size of the "universal graph" which contains all graphs (and 
which therefore has the largest possible orbit size) and arises by somehow renormalizing 
the arising infinite pre-factor. 

The justification for the prefactor in Vrnff = a {[l})VDiff is even more serious given the 
superselection principle derived in fllPf . Let us shortly review how strongly diffeomor- 



phism invariant observables define a superselection rule : Assume that 7(f) and 7 (J) are 
different graphs and that O is a strongly diffeomorphism invariant observable which is 
densely defined on $. There are an at least countable number of diffeomorphisms (p n 
which leave 7(1) invariant but such that v 9 n 1 (7( < ^)) are mutually different graphs for each 
n. It follows that 

O u := < Tj, OTj > aux =< f% n )T 7 , OTj > aux =< T T , Ufo-^OTj > aux 
= <T I ,OU(<p- 1 )T J > aux 

(2.8) 

and so O has the same matrix element Ojj between Tj and an infinite number of mutu- 
ally orthonormal states. Thus, since O was densely defined we conclude that Ou = 0. 
This is the superselection rule : strongly diffeomorphism invariant observables cannot 
map between cylindrical functions defined on different graphs. In particular, they cannot 
map between graphs which are non-diffeomorphic. Therefore, the inner product between 
diffeomorphism invariant states which arise as group averages from different diffeomor- 
phism invariant sectors are completely unrelated and we may therefore choose each of 
them separately. 

Of course, the existence of the superselection rule rests on the assumption that strongly 
diffeomorphism invariant observables are the only ones that givve rise to physical observ- 
ables inquantum gravity. Now, from [l], |[ we know that the Hamiltonian constraint 
operator does not leave the superselected sectors invariant. In particular, the space of 
solutions consists of distributions which are linear combinations of diffeomorphism invari- 
ant states from different superselected sectors. It follows that any non-trivial physical 
observable, that is, one that maps the solution space to itself, and which differs from 
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the identity operator mixes different super-selected sectors. Clearly, such an observable 
cannot be strongly diffeomorphism invariant. Strongly diffeomorphism invariant opera- 
tors are built purely from electric field operators (as they cannot change the graph) while 
weakly diffeomorphism invariant operators will also depend on the connection through 
loops. We conclude that the restriction to strongly diffeomorphism invariant observables 
was too restrictive and that the superselection sectors are therefore spurious. 

Still there remains the issue which value of the constants a([7]) is the correct one. Is 
there is a physical motivation for = 11 

Theorem 2.2 i) The representation of the algebra of strongly diffeomorphism invariant 
obeservables on Ti is completely reducible with TC S ,TC A as invariant subspaces. 
ii) The inner product defined by < r] Diff f,r] Diff g > D iff = (VDifff)(g) forallf,g G $ ; that 
is, a ([7]) = 1, is uniquely selected by requiring that a) strongly diffeomorphism invariant 
self-adjoint operators on 7i are promoted to self-adjoint operators on Hoiff and b) < 
f,9 >aux=< VDifff,VDiff9 >Diff for all f,ge (H s n $) 7 . 

Proof : 
i) 

Notice that each element (p G G(7(/)) permutes the edges and vertices of 7(1) in some 
particular way. This permutation induces a permutation er of the elements of the set 
A ([/]). Thus, o" G S(n([I])) where S(n) denotes the symmetric group of n elements 
(the set of these permutations does not necessarily contain all elements of S(n) but is a 
subgroup thereof). The vectors Tf are therefore invariant under G (7 ([/])) because they 
are completely symmetric under any element of S (n( [/]))■ 

We conclude that for a strongly diffeomorphism invariant operator O the vector OTf 
is again invariant under G{^f{I)) (recall that a strongly diffemorphism invariant operator 
must map functions cylindrical with respect to a graph to functions cylindrical with 
respect to the same graph). Next, notice that the graph symmetry group acts transitively 
and fixpoint freely on the set of vectors {Ti^YlTi by its very definition (every vector 
can be mapped ino any other). It follows that the induced subgroup S'(n([I])) of the 
permutation group on A([/]) contains a cyclic element a of order (that is, one with 

o-n{[i]) — j(j anc j n = n([I]) is the smallest integer with this property). 

To see this just decompose every a G S'(n([I])) into cycles. Take one element which 
contains a cycle of length k and compose it with another element which contains a cycle 
of length / that involves at least one element which is not involved in the first cycle and at 
least one that is involved in the first cycle. Such an element must exist because the graph 
symmetry group acts transitively on A([J]). The result is an element of S'(n([I})) which 
contains a cycle of length k + 1 at least. Proceeding inductively, the assertion follows. 

Given such a cyclic element a G S'(n([I])) corresponding to ip G G(7(/)), let us 
label the states T kjI in such a way that U((p ) k Tj^ = T Ijk+1 . Let us consider a new (non- 

orthonormal) basis of span{T/ fc }^{'' ) 1 given by Aj^ := Tj^ — Ti^+i, k — 1, .., n([J]) — 1. 
That indeed every Tf k is a linear combination of these new vectors follows from the 
identies £f=i T I:l - kT I:k+1 = E?=iP>,z - T I)k+1 ] and Tjj - T Ijk+1 = A I:l + .. + A Ijk . The 
virtue of this new basis is that it allows us to readily verify that 

n([J])-l 

M/ (I :=( E %)% = 0VI=l,.,n([/])-l 

k=0 
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since U((po)^i,i = with Tj )n ([i])+i := 7/^ and because Ei=i Ai,k = 0. We conclude 

that mTl k = for all k = 1, ..,n([I]). 

We will show that a strongly diffeomorphism invariant operator O must leave the cylin- 
drical subspaces of 7i s and 7i A invariant. To that end let be given arbitrary elements / 
of Ti. s R $ and g of 7i A D $. This means that they are finite linear combinations of the 
form / = Y^iLiSiTj. and g = Y,j=i a jTj jk . where kj E {1, ..,n([Jj]) — 1}. We need to 

show that < /, 6g >=< g, Of >= for all / G 7i s H $, g e H A n $. This is equivalent 
to showing that < Tf,OT$ k >=< T$ k ,OTf >= for all J, J and = 1, ..,n([J]) - 1. 
Since a strongly diffeomorphism invariant operator leaves the graph of a spin-network 
state invariant these equalities follow trivially unless 7(f) = 7 (J) so let us assume that 
this is the case. We now use the fact proved above that Tf is invariant under (7(7(1)). 
This implies that, since 7(f) = "/(J), fjjTf = Tf. Since by assumption [fjj, O] = we 
have 

< Tf, OTf k >=< tfjTf, OTf k >=< Tf, Ov.jT A k >= 
by construction of fjj. Similarily 

< T A k , OTf >=< T A k , drfjTf >=< f,jT A k , OTf >= . 

As an aside we notice that G(j(I)) leaves also Ti A invariant because the identity operator 
is obviously strongly diffeomorphism invariant. 

ii) The most general group averaging with respect to adapted bases is given by f]Diff := 
a([y(I)})r]Diff and the associated most general inner product is given by < ^Difffi VDiff9 >Diff'-= 
(r]jr)hff) (g) where a([y(I)}) > depends only on the orbit of the graph. Let us see which 
restrictions are imposed by requiring properties a), b) stated in part ii) of the theorem. 
We begin with b). We have already seen in Theorem ( |2.1| ), i) that the rjoiffTf are or- 
thonormal with respect to the inner product given by a ([7]) = 1 for all [7]. It follows that 
for any f,g 6Wj we have with / = Ej fi T f, 9 = £/ 9i T f 

< V%ff, n% f 9 >Diff= E Ji9j{v% f Tf){TS) = o([ 7 ]) E fi9i = «(W) < /, 9 >«, 

1, j 1 

which fixes a ([7]) = 1. 

We turn to a). Since a strongly diffeomorphism invariant operator O leaves TC A invariant 
the only non-trivial matrix elements involving T A k are given by {j}D%ffTfi^{pTfy = 

since VDiffT A k = and goij f(OT A k ) (Tj t ) = since 7]DiffOT A k = because the averaged 

vector lies in 7i A again. So on H A the operator O is trivially self-adjoint because after 
averaging it is the zero operator. Thus we need to check only whether or not it is true 
that 

< VDiffOf, VDiff9 >Diff=< rjmfff, VDiffOg > Diff (2.9) 

for any f,g& 7i s . However, this is trivial : Since O cannot map between functions which 
are cylindrical with respect to diffeomorphism inequivalent graphs we may assume that 
/, g are in fact cylindrical with respect to the same graph 7. Since both lie in TC S and O 
leaves T~C S invariant all vectors in ( |2.9| ) have the orbit size of 7. Thus, since O commutes 
with all diffeomorphisms and is self-adjoint on 7i aU x 

< VDiffOf iVmffQ >mff = E <U(<p)Of,g> 

¥>e-D([7]) 
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= Yl <OU( V )f,g> 

¥>e-D([7]) 

= <U(v)f,Og> 

=< VDifff, VDiffOg > Diff (2. 10) 

which furnishes the proof. 
□ 

Notice that it was crucial in the proof for the satisfaction of the reality conditions that 
strongly diffeomorphism invariant observables leave the spaces TC s ,7i A invariant. Inter- 
estingly, the Hamiltonian constraint operator defined in [p], |[ [J maps the space H s also 
into itself because it treats all parts of the graph democratically. 

The condition that fixes the value a([7]) = 1 is that inner products between vectors on 
are preserved under group averaging. This implies that the UmffTf = [T/]/^/n([J]) 
are a natural orthonormal basis on Tloiff- What is the meaning of this condition ? It 
means that the group averaging defined is a partial isometry from 7i aU x to Ti-Diff (partial 
because rjniff has a huge kernel even on TL S ). We can interprete this as saying that a 
reality condition selected the inner product because for a transformation to be a partial 
isometry means that it satisfies a (partial) adjointness condition with respect to the inner 
product to be determined. 

3 The dual Constraint Algebra 

Recall that the action of the Hamiltonian H(N) constraint on $ was inherently equipped 
with a huge amount of ambiguity which had to do with the precise position of a particular 
edge to be attached to a graph. In other words, there is a certain infinite-dimensional 
family of prescriptions p and associated Hamiltonian constraints H P (N) such that if / 7 
is a function cylindrical with respect to a graph 7 then there are projections H p (N)f^ = 
H p ^ (N) f\ = X^ e y( 7 ) N{v)H p where Hffi is an operator which attaches certain edges 
in between pairs of edges of 7 incident at the vertex v in a prescription dependent and 
graph dependent way. However, that prescription is covariantly defined, that is, if we map 
the graph with a diffeomorphism (p then the prescription p(y?(7)) is such that the operators 
U(tp)H p ty and t^^ffl are related by a diffeomorphism which in general differs from the 

identity. This has the important consequence that the dual of the operators H P (N), when 
evaluated on diffeomorphism invariant distributions, are independent of p, at least within 
the same diffeomorphism equivalence class of prescriptions p. Thus, with this restriction, 
the dual operator H'(N) defined by H'(N)^(f) := ^(H p (N)f) is completely independent 
of the choice of p, provided that \& e $' is a diffeomorphism invariant distribution. Now, 
the problem is that while \l/ is diffeomorphism invariant, H'(N)^ is not. In fact we have 
the following theorem which demonstrates that the commutator algebra between dual 
diffeomorphism and dual Hamiltonian constraint closes in the way expected from the 
classical Poisson constraint algebra {H (M), V(N)} = H(C^ M ) where V is the classical 
diffeomorphism constraint. 

Theorem 3.1 The quantum constraint algebra between the dual Hamiltonian constraint 
operator H' (N) and the exponentiated diffeomorphism constraint operator U (ip) is anomaly- 
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free on diffeomorphism invariant distributions ^ , that is, 



([H(N), U(<p)])'V = H'{ip*N - N)ty . (3.1) 

Proof : 

Notice that it makes sense to check the constraint algebra only on diffeomorphism invari- 
ant elements of $' because this is where H'(N) is defined. Moreover, since (U((p))'ty = \I/ 
it makes sense to compute the commutator. 

So let fry be a function in $ cylindrical with respect to a graph 7. Then we have by 
definition 

-(([a(N)M<p)))'*)(fJ = ([&{N)M<p))*)(fj) 

= ([H'(N) - U'( V )H'(N)mf,) = *(lH p (N) - frWUfrMJ 
= E {N(vMH^f,) - N^MH^f^)} 

= (H'(N-^N)^)(f,). (3.2) 

The only non-trivial step in this computation has been the last one in which we have 
used the fact mentioned before that the prescription is covariantly defined and that the 
images of the Hamiltonian constraint on diffeomorphic vectors are diffeomorphic so that 
the difference is not seen by 
□ 

We also recall the following theorem from [[10] . 



Theorem 3.2 The commutator algebra of the exponentiated diffeomorphoism constraints 
closes, that is, 

U{ip)U{ip')U- l {<p)U- l W) = U(<p o ip' o ^ o (cp)')- 1 ) . (3.3) 

Proof : 

The proof follows trivially from the fact that U(ip)f 1 = f^) which in particular shows 

that U is a unitary representation of the diffeomorphism group on 7i. 

□ 

Theorems ( p. 1|) , (3.2) mean that the quantum Dirac algebra closes in the appropriate sense 
between Diffeomorphism constraints among themselves on 7i and between Diffeomorphism 
constraints and the dual of the Hamiltonian constraint on H-Diff- 

What remains to be shown is that the algebra closes between two (duals of) Hamiltonian 
constraints. In fact, the following theorem was proved in [0. 

Theorem 3.3 The quantum constraint algebra of the Hamiltonian constraint closes, that 
is 

(([H(M),H(N))y*)(f 7 ) = (3.4) 
for any diffeomorphism invariant \P e <&', any L G $ and any M, N. 
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Formula ( |3.4| ) is unsatisfactory for two reasons : 

1) On the left hand side one would like to see something like [H'(M), H'(N)} rather than 
([H(M), H(N)])'. 

2) On the right hand side one would like to see something like (/ d 3 x(M^ a N ' — M N ta )q ab Vb)' 
rather than just zero. 

The problem with issue 1) is that the product operator H'(M)H'(N) does not exist on 
'H-Difj- The problem with issue 2) is that none of the operators q ab , Vb exist even on 
7i. Moreover, even if the right hand side would be of that kind, how would one tell the 
difference between this operator and the zero operator ? The point is that whenever 
the algebra is anomaly free then the right hand side must annihilate a diffeomorphism 
invariant state. 

In the sequel we will somewhat improve ( |3.2| ) as follows : 

First we emphasize that the prescription p of |fj] was designed such that not only H'(N)ty 
but also the dual of arbitrary products H p (Ni)..H p (N n ), evaluated on diffeomorphism 
invariant states, for arbitrary iVj is independent of p, that is, there exists a well-defined 
operator [H(Ni)..H(N n )]'. This follows immediately from the fact that in |2j the pre- 
scription was covariantly defined, that is, the prescription applied to graphs that differ 
by a diffeomorphism ip yields graphs that differ by some diffeomorphism <p' which does 
not need to equal <p. We stress this point here one more time for the sake of clarifying 
things, however, it is already contained in This fact can be used to define products 
of dual operators indirectly to be the dual of the product and this was actually done 
in Theorem 3.1 and section 3.4 of JJ. (For a different approach that involves a certain 
enlargement of the space of diffeomorphism invariant distributions and therefore enables 



to define products of dual operators directly see p8[ ). Thus, we are able to remedy issue 

1) above. As expected, the dual operator algebra turns out to be completely Abelian as 
one can immediately see from Theorem |3.3| and as pointed out in [|[] . 

Secondly, it turns out that one can quantize the classical expression / d 3 x(M ta N—MN ta )q ab Vb 
in the sense that one can turn it into a densely defined operator 0(M, N) on 7i. We then 
show that 0'(M, N)^f = for any diffeomorphism invariant \I> G This remedies issue 

2) above. 

In conclusion we can write the operator equality (on Tloiff) [H'(M),H'(N)\ = 0(M,N) 
which is as close as we can get to representing faithfully the quantum Dirac algebra in 
the present framework which suffers from the unavoidable fact that H{N) is state depen- 
dency defined (in other words, A l a is not defined, only holonomies along paths are defined 
and the choice of path is the source of the ambiguity) . 
We begin with issue 1). 

Theorem 3.4 Let ^ G $' be a diffeomorphism invariant state and f G $. Then 
^/(H p (Ni)..H p (N n )f) depends only on the diffeomorphism invariant properties of the pre- 
scription p for any Ni, .., N n and any n = 1, 2, ... 

Proof : 

We can assume that / is cylindrical with respect to a graph 7. Also, the ambiguity is 
already completely encoded in the precription dependence of the Euclidean Hamiltonian 
constraint H E because the Lorentzian constraint is just a polynomial of these with oper- 
ator valued coefficients that leave the graph invariant [[|. 

Now, whenever H P E acts on a function cylindrical with respect to a graph 7 then it re- 
duces to a sum of terms H Evee ,, one for each vertex v G ^(7) and one for each pair 
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e, e' G £^(7) of edges of 7 incident at v. The prescription p assigns to each H P E v e e ,, proper 
segments s C e, s' C e' incident at t> and an arc a incident at the endpoints of s, s' and not 
intersecting the graph anywhere else together with a diffeomorphism invariant condition 
about the topology of the routing of a through the edges of 7 incident at v. The "length" 
of s, s' and the "shape" of a is the whole source of the precription dependence. The point 
is that when H E acts again, the new arcs to be attached will never intersect the old arcs 
already attached because at each stage s, s' are always proper segments of e, e'. Consider 
now a particular term 

HK^rA^A - k E^e'J) ( 3 - 5 ) 

in the expansion of \1/ (H E (N\) . .H E (N n ) f) . If it does not vanish then \1/ is a linear combi- 
nation of group averaged cylindrical functions and one of them, ty', depends on the graph 
7' = 7 U ai U .. U a n and it is the only term of \1/ that contributes to (|5T5|) . Here 7 C 7 
is the graph 7 with possibly some of the s,,s^ removed 0, [3]. Now, since we allow for 
smooth diffeomorphisms that leave a graph analytic, there are smooth diffeomorphisms 
involved in the definition of that leave 7 invariant and 7' analytic and that alter the 
length of Sj,^ and the shape of Oj arbitrarily while keeping the topology 0/7' invariant 
(arcs cannot intersect each other or 7 other than in their endpoints but they can slide up 
and down along e, e' and change their shapes). It follows that ( WM does not depend on 
p but only on its diffeomorphism class [p] . 
□ 

Corollary 3.1 The quantitity [H(N 1 )..H(N n )]'^ defined by ^(H p (Ni)..H p (N n )f) is well- 
defined, that is, only depends on the (fixed) diffeomorphism class [p\. 

Definition 3.1 Let B be the vector space of sums of products of operators H P (N) for 
some p and let B' be the set of dual operators defined by B' = {B'; 3B G B 3 [B'^](f) = 
ty(Bf)\/ty G &Diff, f G $} where again &Diff ^ s the diffeomorphism invariant subset 
o/$'. We turn the vector space B' into an associative algebra by defining a product by 
B' X B' 2 := (B2B1)' . Notice that the product is a well-defined map from B' x B' to B' by 
Corollary 

Corollary 3.2 The algebra B' is Abelian. 
Proof : 

Notice that [i?^,^] = ([-62,-81])' and write ^([B 2 , Bi]f) as sums of expressions of the 
form $>(B[H p (M), H p (N))f) where B G B and /' G $. Now use Theorem Q) in the 
argument of the proof of Theorem ( |3.3| ) (or 3.1 in 0) in order to complete the argument. 
(Essentially, one needs to show that one can slide the arcs along e, e' despite the presence 
of the operator B. This can be accomplished by straightforward application of Theorem 

□ 

Next we turn to issue 2) : In the last part of this section we show that one can in- 
deed define an operator corresponding to the right hand side of {H(M), H(N)} whose 
dual annihilates diffeomorphism invariant distributions. 
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Let u a := MN^ a — M iV and consider the triple of vector fields v a with components 
iy a ) b := 6 b as well as the triple of one forms a b with components (<r b ) a := 5 b . Let Xe{ x , y) 
be the characteristic function of a coordinate box with center x and coordinate volume e 3 . 
The volume of the box as measured by q a b is given by V(x, e) := / d 3 yx t (x, y)^Jdet(q)(y). 
Finally, let V a = ti{F a \ ) E b ) be the classical infinitesimal generator of diffeomorphisms. 
With this preparation we have the following classical identity 



0(M,N) := / d 6 x(u a q ab V h )(x) 



lim — / d 3 x(u a e abc e ijk -^=^)(x) x 
^°4e 6 J -/det(g) 



x jd 3 y{{a%e^e ilm ^=){y) d 3 z Xe (x, z) X e(y, z)((v d ) h V h )(z) 

/det(g) J 



lim — — - — / (to A ){x) x 



det(q) 



e l A e r 



x /((«,») A ^^=)(j,) / d 3 z Xe (^^)x e (y^)(WV 6 )(z) 
'det(g) 



det(g) 



'det(g) 



x / d 3 2x e (a;,z)x e (2/,2;)((wa) b V 6 )(z) 



,. e ijk e Um r {A j ,V(x,e)} A {A k ,V(x,e)} 
lim / (a; A — r J (x) x 



e^o 16k 2 J v V(x,e) 

X /" rf 3 2;X e (x,z)Xe(l/,2)((Wa) b V6)(2) 

lim^^ |(a; A ^(i.e)} A {A fc , ^(x,e)})(x) x 
x / (K) A {A 1 , y^e)} A {A m , y/V(y,e)})(y) X 

X / d 3 Z X e(x,z)Xe(y^)(K) fe H)(^) (3.6) 



where in the first step we just have written 

q ab in 

terms of e l a , in the second step we have 
rewritten everything interms of forms, in the third step we have used the key identity 
II 

2«;sgn(det((e J 6 ))K = {Al(x),V(x,e)} , 



in the fourth step we used that V(x, e)/e 3 — > Jdet(q)(x) as e — > and finally in the last 



step we just used {A a (x), yV(x, e)} = {A a (x),V(x,e)}/(2^V(x,e)). 

Expression ( |3.6|) is just one example of a general theorem proved in [20] which says 
that singular prefactors of Jdet(q) can always be absorbed into Poisson brackets which 
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then are subject to quantization as follows : 

Given a function / 7 cylindrical with respect to a graph 7 we triangulate S in adaption 
to 7 as in . We leave the the quantization V a of V a open at this point and will specify 
it only later. For the moment being we just order it to the left hand side and replace 
Poisson brackets by commutators times l/(ih). The connection A gets replaced, just 
as in H, by a holonomy along the edges of a tetrahedron. We use a basis of su{2) 
with the algebra [rj, Tj] = tijkTk so that ti^A) = — A l /2. Each expression of the form 
tr(rj/i s [/i~ 1 , V l l 2 (v , e)])/ 7 , where v is a vertex of 7 and s a segment of an edge starting at 
v, does not depend on s because it must be gauge invariant at the endpoint of s. Since the 
volume operator only acts at vertices of 7 we find in complete analogy with computations 
in the regulated operator 

6 e (M,N)f, = J2 I 'd 3 z Xe (v,z) Xe (v',z)((v a ) b V b )(z)x 

x ]T e n ™[M(v)N(s n (A)) - M(s n (A))N(v)} x 

v(A)=v 

x tr(r^ Sp(A) [/i;; (A) , V 1/2 (v, e)])tr{r k h a ^ A) [h^ A yV^ 2 (v, e)]) x 
x ]T e rst (a a ) c [s c r {A')-v' c ]x 

v(A')=v' 
X tr^A^^ 

(3.7) 

where in abuse of notation s n (A) denotes also the endpoint of the segment s n (A) whose 
starting point is v(A). 

Fix the value of a,v,v' and let now £ b (z) := (t> a ) b x e (t>, 2;)x e (t>', z). We replace Ef by 
Ef(z) = —ifmS/SA^ in V& and order E to the right hand side. Let f be any function 
cylindrical with respect to 7. Also we abbreviate 

Xl(t) := [h e (0,t)nh e (t, l)U B d/d(h e (0,l))AB 
for each edge e of 7. Then we have 

d 3 xC(x)V a (x)f^ 

£ [ d 3 xe(x)F* b (x) f 1 dt6{z,e(t))e b (t)X i e (t)% 



J2 jyte(e(t))Fi b (e(t))e b (t)Xi(t)^ 



-1 



E limVf(e(4))^(efc))[e 6 fe)-e 6 (ViP:(4)/; (3.8) 



71— >00 

ee£(7) fc=l 



for any partition = to < ii < .. < t n — 1 of [0, 1] which we may choose individually for 
each edge e. 
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We now combine ( |3.7|) and (|3.8|). We choose to take first the e — > limit. (This exchange 
of limits is justified by the fact that one can alternatively also keep e ~ 1/n finite and 
and still arrive at the quantum result below, the virtue being that this way of guiding 



the limit also gives back the original function 0(M, N) on the classical side. See ^ 
for details) . The product of charcteristic functions Xe( v , e (tk))Xt( v ' \ e (^fc)) will vanish for 
any k, n if v ^ v' when choosing e sufficiently small. We conclude that only the terms 
with v = v ' survive in ( |3.7| ). Next, since v is a vertex for sufficiently small e, Xe{v, e{tk)) 
vanishes (keeping n fixed) unless k = because we have labelled edges e such that only 
their starting point can be a vertex of 7. Therefore in the limit e — > the sum over k 
in ( p.8|) disappears and the dependence on e disappears when we replace V(v, e) by the 
operator V v as explained in |[22|| . Finally we contract the indices a in £ b (e(0)) = 5 b from 
([TSD and in (a a ) c [s c r (A') - v c ] = [sJf(A') - v A ] and find that (£70 becomes 

e> n (M, jv)/ 7 = -i 16e ;y EE E x 

" K ueV( 7 ) e6B(7),e(0)=u v(A)=u(A')=i> 

x [<(A) - ^(tO^fa) - v b ]Xl[M(v)N(s n (A)) - M(s n (A))N(v)} x 
x trfo/^A)^^ x 

x trfa/i..^" 1 ^). K' /2 ]) tr (^^(A')K~ t (A')' ^ 1/2 DA ( 3 - 9 ) 

where we have used that the operators in the second and third line do not change the 
graph 7 and X\ = Xl(0) is the right invariant vector field on SU (2) evaluated at h e (0, 1). 
Finally, let us choose t\ for each e individually such that e(t\) coincides with the segment 
of e that appears as an edge Sj(A) of one or several tetrahedrons A of the triangulation 
with v(A) = v and let us call this segment e{t\) := s,(A) =: s(e). Then (|3.9|) only 
depends on the triangulation T. The final step is to replace the curvature expression 
[s"(A') - v a ]F ab (v)[s(e) b - v b ] by h a ( Sr (A'),s(e)) ~ 1 su(2) where a(s r (A'),s(e)) is the loop 
Sx o s 2 S3 1 o sj 1 where Si = s r (A'), S4 = s(e) and s 2 , S3 respectively are some arbitrarily 
chosen segments which do not intersect 7 except in the end points of si,S4 such that 
the loop a looks like a little parallelogramm in the coordinate frame under consideration. 
Moreover, the "arc with a corner" s 2 o S3 1 is subject to the same routing conditions as 
the arcs of tetrahedra (see |J for details). 
Now consider the combination 

[<(A) - v a }F* b (v)[s(e) b - v b ]Xl « tx(h a{Sr{A%s{e)) h e & ^(0,1)) - tv(h e (0, l)d W) ) 

and compare it with the infinitesimal variation of f under the deformation of e = s(e) oe' 
into e = a(s r (A'), s(e)) o e = sios 2 o S3 1 o e'. We find (we suppress all arguments different 
from h e ) 

8f!y ■= fj(he) - fj(he) = fj(h a (s r (A'),s(e))he) ~ f'^K) 
= fj(h e + [/ia(s r (A'),s(e)) - ^SU{2)]h e ) ~ fj(h e ) 

= tr([/l a ( Sr (A'), s ( e )) - lsU{2)]h e dhJf!y + o([h a ( Sr (A'),s(e)) ~ lsT/(2)f) 

= [U((p(s r (A'), s(e))) - id„]/ 7 + o([/i aW A V (e)) - W)] 2 ) (3.10) 

where in the last step we have written the variation of f 1 as the change under a corre- 
sponding diffeomorphism (p(s r (A'), s(e)). Clearly, such a diffeomorphism will not only 
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deform e but also the edge of 7 which has s r (A') as a segment. However, one may choose 
<p(s r (A'), s(e)) such that the corresponding change of /' is of higher order. 

Expression ( p. 10 ) shows that (up to higher order terms which vanish in the continuum 



limit) the operator V& can be multiplicatively regulated by multiplying it with a small but 
finite line element (here given by £ a ) and writing it as a small but finite diffeomorphism 
minus the identity operator. This is because the infinitesimal operator does not make any 
sense on Ti as mentioned earlier and since £ a comes from q ab , it happened what was to 
be expected : While q ab , Vb individually do not make sense as operators, the combination 
q ab Vb in fact does make sense. 

Thus, combining ( |3.9|) and (|3.10|) we obtain the desired end result (T reminds us of the 
triangulation dependence) 

6 T (M,N)f y = -i mei f ilm e rst e n Pq J- E X 

^ K veV(-r)eEE(-y),e(0)=vv(A)=v(A')=v 

x [U( V (s r (A'),s(e))) -id H ][M(v)N(s n (A)) - M(s n (A))N(v)] x 
x trfo/^Ajt/i;^^ x 

x tr^A^;^ . (3.11) 

The only point which is important is that that the diffeomorphism operator 
U(ip(s r (A'), s(e))) stands to the left in ( |3.11| ) which implies that when evaluating ( |3.11|) 



on a diffeomorphism invariant distribution the result vanishes. That such an ordering 
is possible while still obtaining a (precription dependent) densely defined operator is nev- 
ertheless a non-trivial result. 

Moreover, Ot(M, N) is clearly just like H(M) prescription dependent but its dual 0'(M, N) 
is not. We may therefore simply set [H'(N), H'(M)} = 6'(M,N) on & Diff . 
Formula ( |3.11| ) deserves several comments : 



1) The immediate question is : Can we do even better than this ? That is, can one 
recover ( fj . 1 1| ) by simply evaluating the commutator [H(M), H(N)]f^ ? Comparing 
( |3.11|) with the expression found in for the commutator shows some similarities 



but also some key differences : 

The commutator, at least for the Euclidean Hamiltonain constraint, involves two 
factors of operators V while (|3.11|) involves four factors of V 1 ^ 2 . Thus, while the 
overall power is the same, the factor ordering is different. Next, the commutators 
involve two loops but (|3.11|) only one. Checking on generic states f 1 reveals that the 



two loops do not in general combine into only one loop. Trial and error reveals that 
one has a chance to combine the two loops into one if and only if the Hamiltonian 
constraint acts also at the vertices it creates. However, then one finds that when 
evaluating on diffeomorphism invariant states, the result does not vanish any longer. 

This is therefore a major difference which seems insurmountable and shows that the 
commutator and (|3.11| ) are simply very different quantizations of the same classical 



object (meaning that, since the classical expression 0(M,N) can be written as a 
Poisson bracket, we may choose to quantize it through the corresponding commu- 
tator). The deeper reason for this difference is the following : when one does the 
classical computation to prove that indeed {H(M), H(N)} = 0(M,N), one has 
to do several rather intricate algebraic manipulations and integrations by part. In 
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particular, the reason why one curvature term disappears is because the Poisson 
bracket picks the partial derivative of A in F. None of these manipulations can 
have a quantum analogue because F, A are not well-defined operators on 7i. It is 
therefore not surprising that the commutator does not give Finally, it is 

reassuring that the combination [M(v)N(s n (A)) — M(s n (A))N(v)] in Q3.11|) does 
appear in the commutator. 

2) That 0(M, N) can be given a meaning at all is astonishing at first because we know 
that the infinitesimal generator of U (if) as well as q ab do not exist on H,. However, 
since both operators combine to a density of weight one, the structure theorem 
proved in |20| applies and it is therefore clear that the integrated expression 0(M, N) 



can be promoted to a densely defined operator on 7i. 

3) Concluding, we see that the quantum Dirac algebra closes in the appropriate sense 
and gives a faithful representation of the classical Dirac algebra in the maximal 
possible way given the choice of representation TC. 



4 moduli, Superselection Sectors and Separability 

In this short section we mainly wish to express a speculation which needs further inves- 
tigation to make it definite. The issue is the following : 

Notice that we have group averaged the states with respect to (subsets of) smooth diffeo- 
morphisms and that graphs were piecewise analytic. Consider a graph with an n— valent 
vertex v. The diffeomorphism group at v reduces to a GL(3, R) matrix and so one has 
at most d 2 = 9 free parameters to change the tangent directions of all n edges incident 
at v. The number of angles 9 (with respect to some background metric) between those 
tangents is given by n(n — l)/2 which for n > 5 exceeds 9 : This is the first example 
of the appearance of the 9 moduli, that is, there is diffeomorphism invariant information 
encoded in the intersections of analytic edges. These moduli take values in a compact 
subset of IH with the cardinality of R. It is these 9 moduli which make H-Diff & non- 
separable Hilbert space ! Other kinds of moduli which are invariant under smooth, that 
is, C°° diffeomorphisms are the the degrees of differentiability of intersections of analytic 
edges, however, these moduli are integer valued and therefore do no spoil separability. 

The question arises whether these moduli have any physical significance or are just 
an artefact of the fact that we are dealing with piecewise analytic graphs and smooth 



analyticity preserving diffeomorphisms. In fact, in |29| it is shown that if one works, for 
instance, with piecewise linear graphs or simply enlarges the diffeomorphism group, then 
the moduli disappear. Another striking feature of the formalism is that geometric oper- 
ators like length, area and volume which one would expect to depend on those moduli 
actually do not depend on them as pointed out in [^0| which is an immediate consequence 
of the fact that their spectrum is discrete. In fact, the only diffeomorphism invariant 
moduli that play a role are orientations of certain tangent directions and so are of the 



discrete type |31|. The speculation is that none of the physical operators can detect the 



continuous moduli. By this we mean the following : 

Notice that the number of moduli parameters is countable : a concrete counting scheme is 
provided by numbering graphs by some integers corresponding to the number of vertices 
together with their valences (this fixes the number of edges) some integers that fix the 
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routing of edges between vertices and some integers k that label C k properties of differ- 
ent iable intersections of analytic edges (we do not need continuous parameters labelling 
the germs of analytic edges because we average with respect to analyticity preserving 
smooth diffeomorphisms rather than analytic diffeomorphisms) . Let 9 := (6*i, 6*2, --) de- 
note a sequence of such a complete set of moduli. Since functions cylindrical with respect 
to graphs labelled with different moduli parameters are orthogonal with respect to 7i 
they continue to be orthogonal with respect to TCoiff by the very definiton of the group 
averaging proposal. Thus, the Hilbert space TCoiff breaks up into a continuous direct 
sum 

and each of the Hilbert spaces 7^5i// corresponds to a fixed sequence 9, all of them being 
naturally isomorphic to each other. Notice that the Hamiltonian constraint does not mix 
these sectors because it attaches arcs to the graph in such a way that the new vertices 
are tri-valent for which there are no continuous moduli. Moreover, TCpiff is separable 
: we have indicated above how to count graphs and otherwise the set of labellings of 
edges and vertices with spins and contractors also is countable. Thus, the corresponding 
spin-network states provide an explicit orthonormal and countable basis. 

G- Supers election Assumption 

The continuous moduli are supers elected, that is, no physical observable maps between 
states of different sectors H^tf. 

Summarizing, if the B-Superselction Assumption turns out to be correct (and the present 
results indicate that this might be the case) then we may restrict ourselves to a separable 
Hilbert space by fixing a value of 9 := # . This means, for instance, that if an operator 
creates new vertices on a graph with valence five or higher then it must create them with 
the correct value of the moduli corresponding to our choice 9_ . 



5 An inner product on the space of physical solutions 

In || a general algorithm was displayed to construct all solutions of the Lorentzian Hamil- 
tonian constraint (and diffeomorphism constraint) of pure quantum general relativity. A 
huge number of solutions have the surprising feature that they are normalizable with re- 
spect to the inner product for solutions of the diffeomorphism constraint. The treatment 
was incomplete in the sense that no inner product for the solutions to all constraints was 
constructed which is capable to render normalizable also solutions which are not normal- 
izable with respect to the inner product to the diffeomorphism constraint. The reason 
for why such an inner product is not entirely straightforward to construct is because the 
favoured version of the Hamiltonian constraint operator constructed in M is not self- 
adjoint (also a symmetric version of the Hamiltonian constraint operator was constructed 
in [H but it involves a rather unnatural technique which makes it less appealing). This fact 



makes the group averaging technique to construct an inner product as given in [10, 24, 23 
inapplicable. 

On the other hand, arguments first raised in suggest that H(M) must not be self- 
adjoint in order that the algebra closes. In fact, the symmetric operator H'(M) +H'(My 
is easily seen to be anomalous. 
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We therefore conclude that we need to deal with the not self-adjoint operator H(M) 
and are asked to invent a group averaging technique for not self- adjoint operators. The 
following quantum mechanical example will hint how that might be done in the case of 
interest. 

5.1 A two-oscillator example 

We consider the Hilbert space 7i := ^(R, dx) ® ^(R, dy) corresponding to a system with 
two degrees of freedom. H plays the same role as &Diff does in quantum general relativity. 
In order to model the Hamiltonian constraint of quantum gravity it will be sufficient for 
our purposes to consider only its Euclidean part H E (M) which is simpler, at least to start 
with (we will generalize it later to the Lorentzian case). The most important property 
of H E (M) is that it creates new edges when applying it to a function cylindrical with 
respect to a graph 7. Thus, we may think of H E (M) as some kind of creation operator in 
harmonic oscillator terminology where the excitation of the harmonic oscillator is replaced 
by an excitation of the gravitational field along the new edge. We therefore study the 
following oversimplified model operator corresponding to H E (M) 

C:=a<-ctt (5.1) 

where c is an arbitrary complex number and a, b respectively are the usual annihilation 
operators associated with the x, y coordinate respectively. Clearly the operator C is far 
from being self-adjoint on 7i and so we cannot use the usual group averaging procedure 
to construct an inner product on the space of its solutions. 

However, this model is simple enough so that all the solutions to C can be constructed 
in closed form, just as in the gravitational case. The complete solution space will suggest 
a natural inner product theron. 

So, let \m,n > be the usual harmonic oscillator eigenfunctions, that is, d\m, n > = 
y/m\m — 1, n >, a^|m, n >= \/m + l\m + 1, n > so that [a, a)} = 1 and similarly for b. In 
fact, a* almost precisely models the operator tr([h a — /i~ 1 ]/i s [/ij 1 , V]) which is the basic 
building block of the Euclidean Hamiltonian constraint. Here a is a a loop bounding a face 
of a tetrahedron A considered in previous sections and s is an edge of A. Namely, does 
not only rise the excitation by one unit (which corresponds to h^ 1 ) but also multiplies by 
a coefficient y/m + 1 which depends on the "spin" m (corresponding to hjyh' 1 , V]). Even 

the power 1/2 of the spin is correct : The volume operator V is bounded by \J\ J 3 | where 

J is the total spin of the graph. Thus it should have eigenvlues of order m 3 / 2 . Now the 
commutator acts roughly like a spin derivative, leaving us with a leading order of m 1 / 2 . 
The span of these oscillator eigenfunctions constitutes a topological vector space $ which 
is dense in 7i and is the precise analogue of (diffeomorphism invariant) spin-network 
states with "spins" m, n. We are interested in distributions \1/ G $' (the dual of $) where 
$ C H C $' such that ^(C\m, n >) = (C"^)(|m, n >) = for all m, n. We are looking for 
^'s of the form ^(f) := / dxdy^lf and then the dual operator can be explicitly computed 
and is given by C = a — cb. That is, it is a complex linear combination of annihilation 
operators, just as in the gravitational case. We can immediately see why C has a huge 
number of solutions ^ which actually lie in $ : Define the level of an oscillator state 
|m, n > to be the number m + n. Then, if / e $ is a finite linear combination of oscillator 
states none of whose level exceeds a given number N, it follows that (C') N+1 f = because 
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(C') N+1 is a linear combination of monimials of the form a k b N+1 ~ k all of which reduce 
the level of an oscillator state by iV + 1 or annihilate it. In other words, C is nilpotent 
on states of a given level N. 

Explicitly, the complete solution space can be constructed as follows : 
The most general state of level N is given by 

N 

\N,{c}>= Y,Cm\m,N -m> (5.2) 

m=0 

and these states (over)span $ as the c m take all complex values. We apply C' and find 

N 

C'\N, {c} >= ^2 c rr \\pm\m — 1, N — m > —c\J N — m\m, N — m — 1 >] 

m=0 

N N-1 

= V] c m \fm\m — 1, N — m > — c m VN — m\m, N — m — 1 >] 

m=l m=0 

N-1 N-1 

= ^2 Cm+iVm + l\m, N — 1 — m > —c ^2 c mV ~N — m\m } N — 1 — m >] 

m=0 m=0 

N-1 

= [c m+ i\/m + 1 — zc m V N — m]\m, N — 1 — m > . (5.3) 

m=0 

Since C' reduces the level precisely by one unit it follows that the condition (C'ty)(\n, m > 
) = is only a condition on those vectors involved in the definition of ^ with the the 
same level. Thus any solution \I/ is a (not necessarily finite) linear combination of vectors 
of the form fl5.2|) . We find the condition 



c m+ iVm + 1 — cc m V N — m = or c m = c' 

and if we wish to normalize the solution then |co| = (1 + \c\ 2 )~ N ^ 2 . Let us simply denote 
the solution of level N by \N >. Of course not every solution need to be normalizable 
since we check (C'^f)(f) only with elements / G $ so that there is no growth restriction 
on the coefficients in ^ = J2n=) ^nW >• 

Another way of looking at what has happened is that we performed a Bogolubov 
transformation on the annihilation operators and require that the physical vacuum is 
annihilated by one of the two new annhilation operators. The states \N > are therefore 
nothing else than the excited states corresponding to the other (independent) annihilation 
operator. 

We now address the question which inner product to impose on the solutions ^ so 
obtained. Usually one expects that a typical solution to the Hamiltonian constraint is 
not normalizable with respect to the inner product on 7i. The intuition behind this 
expectation is experience with typical solutions of self-adjoint quantum constraints that 
can be obtained by group averaging. Let us give an example : 

Consider the same Hilbert space as above but instead of the constraint C consider the 
constraint H = p—q where p, q are the momenta conjugate to x, y respectively. Thus, H is 
a self-adjoint operator on TC and we may follow the group average proposal to construct the 
space of solutions and an inner product thereon. Denote by \p, q > the usual momentum 



N 
m 



I c 



(5.4) 
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generalized eigenstates. Notice that these kinematical states, in contrast to the \m,n >, 
are already not normalizable. Let now \f >G $. Then the group average map is given by 



l/> := v\f>--= I £e**\ f> 

JR lH 

= dpdq —e ttH \p,q><p,q\f>= dpdqS(p,q)\p,q><p,q\f> 
jr 2 Jr Zn Jr 2 

= / dp\p,p X p,p\f >, (5.5) 

JR 



that is, fj = J dp\p,p >< p,p\. Notice that H\p,p >= 0. 
The group average inner product is defined by 



< f,9 >-:=< f\g >=< f,V9 >= J dp< f\p,p X p,p\g > 



(5.6) 



In particular, we have || \f > || 2 = / dp\ < p,p\f > | 2 < oo since || \f > || 2 = 
J dpdq\ < p,q\f > | 2 < oo (under suitable conditions on < p, q\f > which defines $) but 
II 1/ > II 2 = Jdpdq < p,p\f >< p,p\q,q >< q,q\f >= 6{0)\\ \f > \\t = oo blows up. 
This feature is shared by all models whose self-adjoint constraint operator has continuous 
and unbounded spectrum. 

We will now study a second example with a self-adjoint, unbounded constraint operator 
whose spectrum is, however, discrete. Here we will encounter a surprise : all the group 
average solutions are normalizable elements of the Hilbert space. 

Consider again the Hilbert space H as above but now the self-adjoint constraint operator 
is given by H := a^a — tfb, that is, the difference between two harmonic oscillator Hamil- 
tonians. Clearly, the spectrum is entirely discrete and is given by the set of integers. The 
oscillator eigenstates \m,n > are eigenstates of H with eigenvalue m — n and clearly any 
solution is a (not necessarily finite) linear combination of the states \m,m > for arbitrary 
m > 0. This is precisely what the group average map does as well : we have 

l/> := v\f>--= L^e**\f> 



s 1 2n 

= TT—e ttH \m,n><m,n\f>=J2$m,n\m,n><m,n\f> 

m,n-'- 7T ^ n m,n 

= ^|m, m >< m,m\f > (5.7) 



that is, f) = J2 m \ m i m >< m ) m \- The reason why t now only ranges over S 1 is because the 
spectrum of H is discrete in contrast to the above situation and so the group average must 
produce a Kronecker delta rather than a Dirac delta. Notice that the group averaging 
inner product is given by 

< f,9 >~= E < f\ m >< m \9 > (5.8) 

m 

so that in particular 

II 1/ > \\i = E m I < m,m\f > | 2 < Em,nl < m,n\f > | 2 = || |/ > || 2 and 

II 1/ > II 2 = Em,« < m,m\f >< m,m\n,n >< n,n\f >= 6 0fi \\ \ f > \\l ! 

Not only is the norm of an averaged vector with respect to the average inner product 
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smaller than or equal to the norm of the original vector with respect to the original inner 
product but also the average inner product and the original inner product in fact coincide 
on averaged vectors, nothing blows up. The reason for this is that the spectrum is discrete 
so that generalized eigenfunctions are actually elements of the Hilbert space. 

Let us return now to the constraint operator C above and let us set c = 1 for simplicity. 
We are actually able to describe explicitly the complete spectrum of C', namely, it is the 
entire complex plane. Generalized eigenvectors are given by \z,z' >:= \z > ®\z' > where 
\z > are the usual coherent states and the eigenvalue is given by z — z'. Now, while the 
spectrum is not discrete just as in the case of H = p — q, still the generalized eigenvectors 
are proper, normalized (although not orthonormal) elements of the Hilbert space just as 
in the case of H = a) a — ¥b which is in contrast to the case H = p — q. As we will see, it 
is this difference that makes the group averaged solutions again normalizable with respect 
to the original inner product. 

To see this we must first define a group averaging. To begin with notice that all vectors 
of the form \z,z > are annihilated by C . Moreover, recall that coherent states are 
overcomplete but that still an identity operator is given by 

f dzdz rdz'dz' , r ,_ r ,_, , , 

±T-i — I ; — / ; — \z, z >< z, z \ =: / d/i{z, Z) / dj2{Z ,z)\z,z >< z,z I . 



2in J 2in 

Thus, in analogy with H = p — q we try to define an average map as follows : 

dt 



\f> := V\f>:= ff-e*"\f> 

J Z7T 

dt 



d/i(z,z) J dfi(z',z') J ^e itc '\z,z' >< z,z'\f > 
d\x{z,z) \ d^(z' ',z')5(z,z')\z,z' >< z,z'\f > (5.9) 



where S(z) := 5(x) x=z for complex z denotes the analytic continuation of the 5 distribu- 
tion. Interpreting this formal object as S(x)5(y) with z = x + iy results in 

1/ >:= J dfi(z,z)\z,z >< z,z\f > (5.10) 

that is, 17 = / dfi(z, z)\z, z >< z, z\. Another way to justify (|5.10| ) is to interprete S(z) as 
a holomorphic 5 distribution in the sense / dzS(z, zo)f(z) = f(zo) for holomorphic / and 
to drop the remaining dz integral. A reader who does not like the formal derivation (|5.9|) 
may take ( p,10[ ) as a definition which is very natural and motivated by all the examples 



above. In fact, apart from the explicit measure dfi(z,z), this fj is the most general one 
that maps |/ > to a solution. Notice that fj is the precise analogue of the average maps 
/ dp\p,p >< p,p\ and J2 m \ m , m >< m,m\ encountered above, all of them are projectors 
on the solution space and the measures dfi(z, z),dp,l are the diagonal measures that 
survive from the insertion of the 1^ after evaluating the S distribution. 
It turns out that our t] can be written in a yet much simpler form in terms of the vectors 
\N > above. First, by expanding the coherent states in terms of oscillator eigenfunctions, 
we observe that 



oo „m+n 
13 x ^ Z 



\z,z> = e 1 1 \m,n > 

m ,„= vm!n! 
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oo N 



^2 Z N E \ m i N — m > 

n=o m=o Jm\(N — m)\ 



oo N 

= e -\*\ 2 y ± =2 N / 2 \N> (5.11) 

N=0 ViV ! 

where in the last step we have used definition ( |5.4j) of \N >. Formula ( |5.11| ) looks like 
the coherent state \y2z > except that the Fock states \n > are replaced by the solution 
states \N > of C of level N. This confirms our earlier statement that every solution of 
the constraint is a linear combination of solutions of the form \N >. 
Next we explicitly perform the z, z integrals involved in fj and find (the calculation is 
exactly the same as the one for proving overcompleteness of coherent states) 

oo 2( M + N )/ 2 



Y \M >< N\ / dn(z,z)e- 2 W z M z N 



M,N- 



1 1 r 

- Y , \MxN\ dy(z,z)e-W 2 z M z N 



M,N=0 



1 oo 1 
1 oo 

= -y\N><N\. (5.12) 

/ N=0 

Apart from the factor 1/2 expression ( 5.12p is the expected result : the fj operator is a 
projector on the solution space which can be more conveniently written in terms of the 
\N >'s. In the sequel we will disregard the 1/2. 
The group average inner product becomes 

oo 

<f,9>~= E <f,N><N,g> . (5.13) 

N=0 

Notice that the \N > are orthonormal among each other with respect to the original inner 
product. Let fj' be the projector on the subspace of TC orthogonal to the completion of 
the span of the \N > so that 1-^ = fj + fj'. Then, since both fj,fj' are positive operators 
we have 

lil/>ii 2 ~=< Mf ><</,/>= II l/> II 2 

and 

|| |/> || 2 =< f,tffjf >=< f,fjfjf >=< f,fjf >= || |/ > \\l 

and we encounter the same phenomenon as for the constraint H := a) a — b^b ! 



What we learn from these model investigations is the following : 

1) The assumption that averaged vectors are in general not normalizable with respect to 
the original inner product is false. Roughly, it is false for models in which generalized 
eigenvectors are in fact elements of the Hilbert space. Examples include unbounded 
self-adjoint (with discrete spectrum) and non- self-adjoint constraint operators (with 
continuous spectrum) as we saw above. 



30 



2) In the examples above we saw that the following chain of inequalities holds 

iii/>ii* = iii/>ii 2 <ni/>ii 2 (5.14) 

that is, whenever \f > has a norm, so does \f > and moreover we can compute it by 
means of the original or the group average inner product. Thus, the group average 
map can be extended from $ to all of 7i. 

3) However, it is not true that < f,g >=< f, g >^ for f ^ g. Thus, the group average 
inner product is still needed and different from the original inner product. It is 
given by computing the inner products of the projections of general vectors into the 
solution space. The observation that 17 is a certain kind of projector is a general 
observation and can serve as an abstract definition of the average map for a general 
theory whose constraint operator is not self-adjoint. 

4) One certainly can write down solutions to C' which are well-defined distributions on 
$, however, unless they are actually elements of 7i they are not normalizable with 
respect to the average inner product. All distributions which do not have a norm 
with respect to the original inner product have to be discarded. 

5) The importance of coherent states for the solution of the quantum constraint lets us 
expect that coherent states will also play a quite important role in quantum gravity. 
Notice that coherent states are infinite linear combinations of the basic states \N > 
and so they are not cylindrical, however, they are elements of the Hilbert space 
while the states \N > are in fact cylindrical. 

6) We may view the definition 

fi := / dv(s)\s >< s\, 
Js 

where S is a complete set of labels of (generalized) solutions | s > to a (not necessarily 
self-adjoint) dual constraint operator C' and v is a "natural" measure thereon, as 
a more general definiton of a group average map which in the case of self-adjoint C 
reduces to the usual definition. The measure u, as we have seen, arises typically as 
follows : one considers first the spectral measure \x for the constraint operator C 
and then dv = J s dfi8(C') is the measure induced by \x by deleting the integral over 
the spectrum S of C. In case of a self-adjoint constraint operator this is precisely 
what happens via the traditional group average approach. 

Finally, let us address the question of how to model the Lorentzian Hamiltonian constraint. 
Recall from []2| that the Lorentzian Hamiltonian constraint consists of two pieces H(N) = 
T(N) — H E (N) where H E (N) is the Euclidean Hamiltonian constraint which can be 
modeled by a linear combination — cb as above. The piece T(N) on the other hand is a 
complicated multiple commutator between two factors of H E (N) and three factors of the 
volume operator V and several holonomy operators along open paths. More precisely, it 
contains two factors of the form hj 1 [h s , [V, -£^(1)]] and one factor of the form h~ 1 [h s , V] 
where s is an open path. The operator V does not change the graph, it acts essentially 
by multiplication with functions that depend on the spin. The above combinations of 
holonomies along open paths and back also do not change the graph. Thus it is only the 
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two factors of H E (N) that change the graph. Moreover, since every commutator decreases 
the spin power by one unit we conclude that the spin power of h~ 1 [h s , [V, H E {1)]] is of 
order zero (recall that the spin power of H E (N) is 1/2 to see this). Thus the spin power 
of T(N) is 1/2 as well. Since the "spin" power of a) is 1/2 as well we conclude that 
the bilinear operators (fi^) 2 , (b^) 2 , appropriately model T(N) : they create two new 
excitations and their spin power is even higher than that of T(N). In order to account 
for the difference in the spin power one could choose to divide by a square root of the 
number operator but we will not do that because that wouuld pushh the analogy too far 
and even if we would do it the argument displayed below would still apply. We therefore 
propose to study the constraint operator 

C = a(a 2 + b 2 ) + f3ab - 7(0 + 6) 

where the choice of coefficients respects the fact that the coordinates x, y have to be 
treated on equal footing in order to mirror the expression of T(N). 

We can immediately solve C : consider a state \z,z' > as above, then C'\z,z' >= 
[a(z 2 + (z') 2 ) + fizz' — j(z + z')]\z, z' > is an eigenfunction of C', again the operator C' 
is diagonalized by proper elements of TC ! The condition that the eigenvalue vanishes 
leads to a condition z 1 = h(z) where h(z) is a holomorphic function of z (possibly on 
a two-sheeted non-compact Riemann surface depending on the choice of the parameters 
a, fi, 7). The same arguments as above now suggest to define the average operator by 

fj :— J dfi(z, z)\z, h(z) >< z, h(z)\ (5.15) 

Expression (|5.15|) is sufficient to prove that for any cylindrical state \f > (that is, they 
are finite linear combinations of the states \m,n >) the averaged state \ f >:= rj\f > is 
still an element of the Hilbert space. Namely, since || \z,h(z) > || = 1 we have by the 
Schwarz inequality 

lll/>H 2 

< J dfi(z, z)\ < f\z, h(z) > I J dji{z' , z')\ < f\z, h(z') > \ \ < z, h(z)\z' ', h(z') > | 

< (fd^(z,z)\<f\z,h(z)>\) 2 . (5.16) 



But < f\z,h(z) >= p{z,h(z))e-^ z \ 2+ \ h ^ where p(z,h(z)) is a polynomial depending on 
|/ > so that ( |5.16| ) obviously converges. We also have 

< f,9 >-=< f,fjg > and || |/> \\t = Jdfi(z,z)\ < f\z,h(z) > \ 2 < 00 (5.17) 

for the same reason. This time we cannot find an obvious relation between || \f > 
IUII 1/ > IN II 1/ > II though. 

We conclude that the normalizability of averaged vectors with respect to the original inner 
product carries over to the more complicated "Lorentzian" constraint. 

5.2 A physical inner product for quantum general relativity 

Let us first check whether the inner product that results from averaging the diffeomor- 
phism constraint actually is of the projector type mentioned above. This is easily seen to 
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be the case : if we write r) Di ff := J2[Tf][Ti][Tj]^ then we recover indeed f) Di ff(Tf) = [Tf], 
i.e. we have demonstrated that T]Diff is of the general type encountered in the examples, 
namely it is a projector formed by elementary solutions which themselves are averages of 
a suitable basis of the original Hilbert space. 

Next, let us recall the structure of the solution space of Lorentzian quantum gravity 
from ||. Consider all possible labels / of spin-network states and call the resulting set 
W. The set W splits naturally into two parts : a piece W , called the set of sources and 
its complement Wo in W. The set Wo is characterized by the property that no vector 
7/ for / e Wq can occur in the (decomposition into spin-networks of the) image of the 
Euclidean Hamiltonian constraint H E (N) on $ for no choice of prescription p. The set 
Wo therefore contains the labels which are sufficient to span the image of H E (N). It can 
be further decomposed as follows : for each I e Wo let W^ (I) be the set of labels that 
one finds by decomposing [H E (N)] n Ti into spin-network states for all possible N and for 
all possible prescriptions p. These are the states of level n with source I and they are the 
precise analogue of the states \m, N — m > of level N considered above. We now take the 
group average with respect to the diffeomorphism constraint and it then turns out that 
all the states with labels [J] coming from different W^ n '(I) (up to a diffeomorphism) are 
orthonormal. Thus the diffeomorphism invariant Hilbert space splits as 

timff = ®[i],iew a T~t[i] where H\i\ = ©£U^[w(")(/)] 

and where [W' n )(I)] denotes the group averaged labels of W^ n \l). The space T~L\w( n )(i)] is 
obviously a finite dimensional vector space and so TC[i] is separable. Moreover, H'(N) maps 
7i[r\ into itself, more precisely, it maps Wrcy(n)rni into Wrcy(n-i)mi UHmrMjmi. Therefore 
the Hamiltonian constraint can be solved on each Tim separately. Let us label the vectors 
of 7i\wW(i)] by (notice that all these vectors are diffeomorphism invariant). Then there 
are matrices afj, bfj which depend on [/] and on one of the (set of diffeomorphic) vertices v 
of the graph underlying I such that H' V T™ = a^T™' 1 — b^TJ 1 ' 2 . Writing a general element 
of $' as a linear combination of distributions of the form (one for each [/]) ^ n d^T^ we 
find that c™ +1 a™ +1 = c™ +2 6™ +2 for each [I],n,j, v. It turns out that the matrices a n , b n are 
sufficiently degenerate as to allow for solutions c™ such that c" +1 a^ :+1 = dl^b™^ 2 = for 
each [I] , n, j, v and so does not mix levels but a general solution will mix levels. 
The precise growth of the coefficients of the matrices o^pb^ is unknown (although the 
reasoning from above suggests that they do not grow worse than j 1//2 where j is the total 
spin of a state). On the other hand, if b was identically zero as it is the case for the 
Euclidean Hamiltonian constraint then levels do not get mixed and then the example 
C' = a — b considered above suggests (although we have no proof at present, of course) 
that a general solution will be normalizable with respect to the inner product on Ti-Diff- 
The example C' = a{a 2 + b 2 + (3ab — 7(6 + b) studied above suggests that the same is true 
for the full Lorentzian constraint, i.e. b 7^ (again we have no proof for this). 
Suppose now that we have explicitly solved the infinite number of linear equations c" +1 a™ +1 = 
c" +2 6^ f 2 for each [I] , n, j, v exactly (this is actually possible, although quite complicated; 
we expect that one can invent some efficient code that lets us write down the solutions at 
least up to a controllable error). The examples studied above suggest that every solution 
is a kind of coherent state. The first possibility is that every solution lies mTi,Diff- Then, 
since H.\i\ is separable, one can in principle find an orthonormal basis {T[j]^}^L (with 
respect to < ., . >z>iff) on the solution space for each [/] which defines actually a subspace 
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of If not every solution lies in TH-Diff then we can still find a countable generalized 
basis {TJ/j^I^Lq built from infinite linear combinations of elements of H.[i\ which are or- 
thonormal (with respect to < ., . >Diff) in the sense of generalized eigenvectors, that is, 
< Tjj] i£ , >Diff— unless [I] = [J] and fj, — u in which case this quantity diverges. 
(This is quite similar to the momentum eigenfunction "normalization" in the sense of 5 
distributions : < p,p' >= 5(p,p'). It is also precisely what happens with the group aver- 
aged spin-network states : with respect to the inner product on 7i they are orthogonal in 
the sense that 

<[T / ],[T J ]>:= J2 [Ti}(T) = 

re{Tj} 

unless [I] = [J] in which case this quantity blows up). 

In either case the group average proposal then leads us to define the projectors 

oo 

VHam ■= W] where t) [7] = T [i],» T \i},^ ( 5 - 18 ) 

[I] fM=0 

physical states as arising from group averaging are defined by T p h ys := fjT for each T £ 
&Diff an d the physical inner product between such physical states is given by 

< T phyS} T' phys > p hys'-=< T, f)T' >Diff = T phys (T') (5.19) 

where the latter notation means evaluation of the distribution T p f iys £ i^'oiff)' on the 
"test function" T' £ ®' Di ff by means of the inner product of TCniff- 

Equation ( fTT9|) is a physical inner product for Lorentzian quantum general relativity 



which is naturally suggested to us by group average kind of reasonings. It is precisely of 
the structure of the group average map for the diffeomorphism constraint rj^iff displayed 
at the beginning of this section. 
Remarks : 



1) We want to stress that the inner product ( |5.19|) exists even if the solutions Tjj] i(U 



are not normalizable with respect to < ., . >Diff although we expect this to be 
true. We have struggled to list arguments in favour of normalizabilty of Tfn j(U only 
in order to remove the feeling of discomfort : The quantum mechanics example 
H = a) a — b*b shows that physical states as obtained by rigorous group averaging 
can be normalizable with respect to the original inner product without that anything 
is wrong with that. 

2) The fact that solutions to the Hamiltonian constraint are presumably normalizable 
with respect to the diffeomorphism invariant inner product can be interpreted as 
saying that the whole distributionality of a solution to all constraints is already 
captured by averaging with respect to the diffeomorphism constraint alone. In fact, 
if we had not split the solution of diffeomorphism and Hamiltonian constraint into 
two steps but had solved all constraints in one stroke without even constructing 
fjDiff then we would not have even noticed this fact. 

3) As one can show, the "flat state" 6(F) can be given rigorous meaning as a distribu- 
tion on &£)iff IHH an d it solves H E (N) at least. In 2+1 Euclidean gravity the flat 
state is a physical state for the curvature constraint F — 0. On the other hand, to 
make 2+1 Euclidean gravity resemble 3+1 Lorentzian gravity one can replace F = 
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by Diffeomorphism and Hamiltonian constraints One finds that in 2+1 it is in- 
deed true that physical solutions as sleeted by the so defined Hamiltonian constraint 
(analogously constructed as in this paper for 3+1 Lorentzian gravity) are normal- 
izable with respect to < ., . >Diff- It follows then that the flat state 5(F) is not 
normalizable with respect to < ., . > p h ys because it is an uncountably infinite linear 
combination of diffeomorphism invariant spin-network states with all coeeficients of 
order unity [^| and so is not even normalizable with respect to < ., . >mff- Thus, 
the inner products of a topological quantum field theory (like 2+1 Euclidean gravity 
in the Witten formulation) and quantum general relativity in 3+1 dimensions have 
inner products which are presumably very singular with respect to each other as 
one would expect. To prevent wrong conclusions, notice that F = 0, which can be 
solved, for instance, by A = 0, does not mean that we are in Minkowski space : 
A = T + K = just does not have any obvious Minkowski space interpretation. 

4) Notice that since the projector inner product ( |5.19| ) is induced from that of 7i for 
which the classical reality conditions are implemented as adjointness relations, we 
can claim that (|5.19p incorporates the physical reality conditions on Dirac observ- 
ables provided a Dirac observable is self-adjoint on 7i already. This is because a 
Dirac observable leaves the solution space invariant and is therefore already pro- 
jected ("it commutes with the group averaging"), see below. 

5) A complete set of Dirac observables in the sense that it be a self-adjoint and densely 
defined operator on TC P h ys are trivial construct : Define 

O'j^ ■= and 6 [IMJ]U := ± (5.20) 

where the dagger operation is with respect to < ., . > p h ys then linear combinations O 
of symmetric operators of the form ( |5.20|) form a complete set of Dirac observables 



(modulo domain questions). 

Notice that in case that the Tm^ are not < ., . >Diff normalizable then ( |5.20| ) has 
to be understood in the folowing sense : there are pre-images G $ such that 
Vf[i]fj. = Now the T[/] M become normalizable with respect to < ., . > p h ys be- 

cause by definition < TW, Tjj]„ > p h ys '-= T^(f^ u ). In case of general relativity the 
are related to the "source states" constructed in || which are really elements of 
$ and orthonormal with respect to < ., . > (after diffeomorphism group averaging 
they form an orthonormal system (but not a basis) contained in Tim)- This implies 
that the T[j] M are orthonormal with respect to < ., . > p h ys - 

From this follows then immediately that r)0 = Of). Also, O is self-adjoint on TC if 
we define < f,0' [IMJ]l/ g >:= T m Xf)T[j\v(g) ■ 
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